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We study the effect of restart, and retry, on the mean completion time of a generic process. 
The need to do so arises in various branches of the sciences and we show that it can naturally be 
addressed by taking advantage of the classical reaction scheme of Michaelis & Menten. Stopping 
a process in its midst—only to start it all over again—may prolong, leave unchanged, or even 
shorten the time taken for its completion. Here we are interested in the optimal restart problem, 
i.e., in finding a restart rate which brings the mean completion time of a process to a minimum. 
We derive the governing equation for this problem and show that it is exactly solvable in cases of 
particular interest. We then continue to discover regimes at which solutions to the problem take on 
universal, details independent, forms which further give rise to optimal scaling laws. The formalism 
we develop, and the results obtained, can be utilized when optimizing stochastic search processes 
and randomized computer algorithms. An immediate connection with kinetic proofreading is also 
noted and discussed. 

When engaged in a specific task for a time period 
that extends beyond our initial expectations, we are con¬ 
stantly faced with two alternatives—either keep on go¬ 
ing or stop everything and start anew. Every now and 
then we opt for the latter, hoping that a fresh start will 
break-off an unproductive course of action and expedite 
the completion of the task at hand. This decision could, 
however, turn out to be counter-productive—nipping an 
awaited, but unforeseen, finale in the bud. To restart, or 
not to restart, that is therefore the question. 

Not at all unique to our everyday lives, a “dilemma” 
similar to the one described above is relevant to virtu¬ 
ally any physical, chemical, or biological process that can 
be restarted. Most notably, restart (or unbinding) is an 
integral part of the renown Michaelis-Menten Reaction 
Scheme (MMRS) illustrated in Fig. 1 [1]. Originally 
devised to describe enzymatic catalysis, the MMRS has 
attracted on-growing scientific interest for more than a 
century [2]. Indeed, nature is full with an astonishing va¬ 
riety of Michaelian processes. DNA-DNA hybridization, 
antigen-antibody binding, and various other molecular 
processes can all be described by the MMRS [3]. That 
and more, the simplicity and generality of the scheme 
have rendered it widely applicable and it is now used to 
describe anything from heterogeneous catalysis [4-6] to in 
vivo target search kinetics [7]. As a matter of fact, one 
can easily convince himself that any first passage time 
(FPT) process [8]—be it the time to target of a simple 
Brownian particle, or that of more sophisticated stochas¬ 
tic processes [9-12] and random searchers [13-15]—can 
become subject to restart [16-22] and is then naturally 
accommodated by the MMRS. Wishing to acquire a uni¬ 
fied view on restart phenomena we identify the MMRS 
as an ideal object of study. 

Central to our understanding of the MMRS is the 



Figure 1: Color Online. The Michaelis-Menten reac¬ 
tion scheme. An enzyme E can reversibly bind a substrate 
S to form a complex ES. The substrate can then be converted 
by the enzyme to form a product P or, alternatively, unbind. 
The conversion of the substrate to a product is manifested via 
the process of enzymatic catalysis. Following either catalysis, 
or unbinding, the enzyme is free to act on additional substrate 
molecules. 


Michaelis-Menten (MM) equation [1[. This equation pro¬ 
vides a hundred-year-old prediction by which any in¬ 
crease in the rate of unbinding (restart) will inevitably 
slow down the rate of enzymatic turnover or, equiva¬ 
lently, prolong the completion time of any process that 
falls into the MMRS category. Surprisingly, this predic¬ 
tion was never tested experimentally, but rapid advance¬ 
ments in single-molecule techniques [23-25] have recently 
motivated us to question it from a theoretical perspec¬ 
tive [26]. Contrary to the classical result, we have found 
that unbinding may also facilitate the successful comple¬ 
tion of a reaction. In the emerging picture, a non-trivial 
solution to the restart dilemma is given by an optimal 
unbinding rate which strikes the right balance between 
the need to abort prolonged reaction cycles and the need 
to avoid premature termination of ongoing ones. Similar 
observations were made in the context of search processes 
[17] and, in particular, in the context of DNA search [16] 
where the authors analyzed facilitated diffusion [27] from 
a very general perspective. The optimal restart rate de¬ 
pends, however, on the full distribution of the underlying 
FPT process (catalysis in the case of enzymatic reactions) 
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[26], and the question of what can be said about it in the 
general case remained open. 

In this letter, we address the optimal restart problem 
within the framework of the MMRS. First, the govern¬ 
ing equation for this problem, Eq. (3) below, is derived 
and solved exactly in several cases that are of particular 
interest. We then show that, in the general case, there 
are two regimes at which solutions to the problem are 
universal. These solutions are given in Eqs. (4) and (5). 
The applicability of our approach is widespread as it al¬ 
lows one to incorporate restart into an existing, generic, 
FPT problem in an almost plug & play manner. 

The optimal restart problem. In formulating the 
optimal restart problem we adopt the terminology of en¬ 
zymatic reactions (Fig. 1) and consider a scenario in 
which the processes of binding, unbinding and cataly¬ 
sis are all stochastic [28]. This probabilistic view point, 
whose origins can be traced back to the work of Ninio [29], 
has found one of its prime applications in the analysis of 
single molecule experiments [23, 30-32], and is now well 
established theoretically [33-35]. Defining the turnover 
rate ktum as the reciprocal of the turnover time—the 
mean FPT required to complete the reaction cycle—we 
will be interested in the turnover-unbinding interplay. 

When binding, unbinding, and catalysis times are 
exponentially distributed with rates kon[S] ([S'] being 
the concentration of the substrate), koff, and kcat 
respectively—the single molecule MM equation is at¬ 
tained [33] 


k 


turn 


kcatiS] 

[S] + Km ’ 


( 1 ) 


with Km = (koff + kcat)/kon (note that, in contrast 
to other rates in the MM equation, the turnover “rate” 
is not a rate in the “exponential sense”). In this case, 
the memory-less property of the exponential distribution 
asserts that the time remaining till the completion of an 
ongoing catalytic step, given its age, is exponential and 
statistically identical to that of a newly started catalytic 
step. It is therefore clear (see Km above) that ktum is a 
monotonically decreasing function of the rate koff- 
Non-exponential time distributions are, however, quite 
common in a variety of complex systems [36-40], and it 
has recently been recognized that enzymes are no excep¬ 
tion in that regard [23, 32, 41, 42[. This result is per¬ 
haps not surprising as catalysis is intrinsically coupled to 
the enzyme’s internal degrees of freedom via a complex 
energy landscape [43] which can give rise to strong devi¬ 
ations from exponentiality and other anomalies [44-48]. 
Renewal theory can then be invoked to provide a gen¬ 
eralized mathematical treatment of the MMRS [26]. A 
completely general analysis of the turnover-unbinding in¬ 
terplay is then very hard, but progress can be made if one 
narrows down to the case of exponentially distributed un¬ 
binding times [26, 33]. Letting koff denote the unbinding 
rate (assumed to be independent of the catalytic process). 


and fcat{t) the probability density function (PDF) of a 
generally distributed catalysis time Teat, it is possible to 
show that [26, 33] 
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where fcat{k) = f is the Laplace transform 

0 

of fcat{t), and {Ton) is the mean of a generally distributed 
binding time. Equation (2) extends the classical result 
of Michaelis & Menten and brings new physics. Indeed, 
an interesting corollary of Eq. (2) is the possibility of 
restart-facilitated-turnover, i.e., a regime in which un¬ 
productive unbinding events lead to accelerated turnover 
[26]. This type of counter-intuitive behavior is categor¬ 
ically precluded by the classical MM equation and is 
therefore considered “non-classical” or “anomalous”. 

In Fig. 2A we use Eq. (2) to plot ktum as a function 
of koff for different catalysis time distributions (CTDs). 
An asymptotic decay of ktum to zero at large koff di¬ 
rectly follows from Eq. (2), is common to all plots, and 
is therefore not shown. At intermediate koff, however, 
unbinding can be either inhibitory {dktum/dkof f < 0) 
or excitatory {dktum/dkof f > 0) and the surprise comes 
from the fact that the latter implies the breaking of the 
classical limit for maximal turnover rates: ktum (0) = 
{{Ton) + {Teat)) ^■ Maximal turnover rates, and the un¬ 
binding rates at which they are attained, can however 
vary considerably. What therefore determines if unbind¬ 
ing will enhance turnover, and what sets the maximizing 
unbinding rate 

The fundamental equation of optimal restart. In 

order to address the questions presented above we first 
derive a governing equation for the optimal restart prob¬ 
lem. Namely, we show [49] that the stationary points of 
the turnover rate are the solutions of 


_ fcat{k){feat{k) 1 ) ^ 1 ^ 
k?dfeat{k)/dk k 


{Ton) , 


( 3 ) 


where the function 'l'(fc) is uniquely determined by the 
CTD. Moreover, we show that 4>(A:) has the following 
property: dktum/dkoff > 0 'l'(fc) > {Ton) and vice 

versa. Consequently, a local maximum of the turnover 
rate is attained at an unbinding rate, , which satis¬ 
fies -^{k/^ff) = {Ton) and -^'{k/^ff) < 0. 

Equation (3) clarifies the role of binding in our prob¬ 
lem. Restart and the initiation of a new turnover at¬ 
tempt, inevitably involve a “penalty”—the necessity to 
go through the binding process all over again. Taking 
the perspective of turnover rate maximization, lengthy 
binding times are hence a deterrent against restart and 
rapid binding an incentive to it. In perfect accord with 
this intuition, we note that a maximum in ktum will de¬ 
velop if and only if {Ton) drops below the critical value 
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Figure 2: Color Online. A. Turn-over rate vs. unbinding (restart) rate for different catalysis time distributions 

[49]. Contrary to classical theory, non-monotonic dependencies are possible and the classical upper limit on turnover rates 
(dashed) can be broken. The distributions are numbered by: 1-Log normal, 2-Double exponential, 3-Weibul, 4-Double Erlang, 
5-Double exponential, 6-Exponential. Heterogeneity is observed despite the fact that {Ton) = 0.1 in all instances, all underlying 
time distributions share a mean of {Teat) = 1, and some even share the same variance {a^{Tcat) = 2 for distributions No. 
1-3). Maxima in turnover rates are denoted by full circles. B. A graphical illustration of the fundamental equation of 
optimal restart (Eq. 3). Extrema of the turnover rate are obtained at points where ^'(fco//) intersects the mean binding 
time (Ton). Maxima from panel (A) are once again denoted by full circles. Lines for distributions No. 5-6, where 4'(fc) < 0, 
are not drawn. Note that line No. 4 intersects {Ton} twice, first with a positive slope (minimum of ktum) and then with a 
negative slope (maximum of ktum)- Inset B. vs. mao:/{Ton)- When {Ton) > k^/f = 0. As {Ton) drops 'I’max, 

a maxima may develop gradually (lines No. 1-3) or abruptly (line No. 4). As {Ton} approaches zero, k^ax can either diverge 
(lines No. 2 and 3) or asymptotically converge to a plateau (lines No. 1 and 4). 


of '^max = rnax{'i>{k)}. Indeed, since ktum is posi- 
fc >0 

live and asymptotically decays to zero as koff —f oo, 
lE'(fc) will eventually intersect any level in the range 
'^max > (Ton) > 0 with a negative slope (provided 
'^max > 0) [50]. In particular, a maximum in ktum 
will develop whenever {Ton) drops below T'o = ^'(0)— 
an observation that will come in handy later on. These 
properties of 'l'(fc) are graphically illustrated in Fig. 2B. 

A classical example and a family of exactly solv¬ 
able cases. The importance of Eq. (3) cannot be 
overstated as it allows one to find optimal restart rates 
for generic FPT processes. In order to demonstrate the 
power of this formalism, we will now reanalyze a prob¬ 
lem studied by Evans & Majumdar in [17]. Consider a 
particle searching for a stationary target via one dimen¬ 
sional diffusion. Setting the initial distance between the 
particle and target to L and the diffusion coefficient to 
D, it has long been known that the mean FPT to the 
target diverges. What happens, however, if the particle 
is returned (restarted) to its initial position with some 
given rate koff (assume (Ton) = 0)? Well, since the FPT 
distribution of the original problem is known to be given 
by fcat{k) = (Laplace space representation of 

the Levy-Smirnov distribution), it immediately follows 
that the mean FPT of the restarted problem is given by 
in Eq. (2)—and we further note that it is finite 
for any positive restart rate! In fact, this is true for any 
fcatik), and as long as {Ton) is finite, but is particularly 


striking when the underlying FPT process is equipped 
with an infinite mean. Moreover, by solving Eq. (3) one 
can see that = {z*)"^D/, where z* ~ 1.59362... is 
the solution to z 12 = 1 — e~^. Clearly, the same modus 
operand! can also be used to study the effect of restart 
on many other FPT classics [8]. In particular, one could 
readily generalize the above example for the one sided 
Levy distribution fcat{k) = (0 < a < 1) to 

obtain k))fY = where z* is the solution to 

az = \ — . 

Analytical solutions to Eq. (3) are hard to find. 
It is thus interesting to note that whenever the 
Laplace transform of the CTD has the following form: 
fcat{k) = (1 -I- afc)/(l + hk + ck"^) (for some con¬ 
stants a, b, and c)—Eq. (3) reduces to a quadratic 
and is hence exactly solvable. A particular exam¬ 
ple in this category is the Exponential distribution, 
fcatit) = (A > 0), for which a = c = 0, 6 = 

and 'l'(fc) = 0. As another example, think of the 
Double-Exponential distribution, fcatit) = pXie~^^*^ + 
(1 — p)X 2 e~^^* (0 < p < 1, Ai > 0, A 2 > 0), for which 

«= - 1 ^). &= + c= and'I'(fc) = 

f P>-l+('t—p)X2 7,2 _|_ 2(AiA2) 7 I AiA 2((1—p)Ai+pA2) \ 

(T-pMTWMP (1-p)p(Ai-A2)2 {1-p)p(Ai-A2)2 J 

Finally, and perhaps most importantly, consider the 
class of distributions which do not fall into the above- 
mentioned form, but can rather be asymptotically 
approximated by it. As we hereby show, the basin of 
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attraction for this class is wide—rendering asymptotic 
solutions to the optimal restart problem (almost) 
universal. 

Universal behavior at {Ton) ~ 'I'o- When {Ton) 
approaches 'I'o, “small /c” solutions to Eq. (3) are antic¬ 
ipated provided 'I'o > 0 (see Fig. 2B). One can then try 
and approximate 4'(fc), at small k, considering that in 

this limit foat{k)^lT E , where M„ = (T"^^) 

n—1 

is the n — th moment of the CTD. It can then be verified, 
utilizing the definition of 'I'(A:), that this expansion must 
be carried out to third order in k (m = 3)—if it were 
to correctly capture 4'(fc) to first order. However, under 
direct substitution of such an expansion into 'l'(fc), Eq. 
(3) becomes a forth order equation and further analytical 
advancement becomes extremely cumbersome. 

To circumvent this difficultly, we make use of the 
widely applied Fade approximation scheme [52] and try 
fcat{k) ~ (1 -I- afc)/(l + bk + cfc^). Doing so, we note 
that this approximation: (i) can be made exact to third 
order in k by proper choice of the constants a, b, and 
c [49]; (ii) decays to zero as /c —>■ oo—as required from 
a Laplace transform, and in sharp contrast to the di¬ 
vergences of any power series expansion; and (iii) ren¬ 
ders the solution to Eq. (3) immediate as it gives for 
k « 1 [49]: 'I'(A:) cs , where 

^<0 = (Ms - 2M2)/2 Mi and Rq = ■ We now 

see that in this limit solutions to the optimal unbinding 
problem are insensitive to fine details of the CTD as they 
are governed by 'I'o and Rq only. 

As we have previously observed, the introduction of 
unbinding is asserted to speed up turnover whenever 
{Ton) < ^0 since this implies 3^"™ > 0. The 

newly derived expression for 'I'o allows us to interpret 
this result probabilistically. Indeed, setting koff to zero, 
it is easy to see that {Ton) < 'I'o, if only if, the mean 
duration, {Ton) + {Teat)-, of a new turnover cycle drops 
below the mean residual duration, f (Toat) / {Teat), of an 
ongoing catalytic step [53]. Unbinding will then have an 
excitatory effect but two distinct scenarios should never¬ 
theless be told apart. 

When Rq > 0, 4''(0) < 0 (e.g., lines No. 1-3 in Fig. 
2B) and, as {Ton) approaches 'I'o from above and crosses 
over to its other side, a maximum of the turnover rate 
gradually develops at 


In particular, setting A = {Ton)~^ — 'f'o we observe 
that to first order ~ 2^- This characteristic de¬ 
pendence is further discussed in Fig. SI [49]. 

On the other hand, when i?o < 0 (e.g. line No. 4 
in Fig. 2B), 'I''(0) > 0, and 'l'(fc) has a local maxima 


f'eatiO) < -/?at(0) 


f'eatiO) > -/Lt( 0 ) 



-1 Diverges to Infinity 0 Starts From Zero 

Starts From 
Positive Constant 


fcatW-t" 
(t« 1) 


Figure 3: Color Online. Asymptotics of the optimal 
restart problem, at fast binding times, is governed by 
the behavior of the catalysis time distribution near 
the origin. 


at some k* > 0. Then, as {Ton) first hits 'I'(fc*) from 
above, both a minimum and a maximum of ktum abruptly 
appear (see “jump” in k)^^^, inset of Fig. 2B). As {Ton) 
continues to decrease, these two extrema drift apart and 
it is important to observe that the small k solution to Eq. 
(3) is then a minimum, rather than a maximum. As {Ton) 
drops below this minimum necessarily disappears— 
leaving behind a maximum of ktum at a point which 
is strictly separated from zero. Before moving forward, 
we note in passing that i?o < 0 if and only if the residual 
duration of an ongoing catalytic step has a coefficient of 
variation that is smaller than unity. 

We end this section by noting that when the catal¬ 
ysis time distribution is “heavy tailed”—as happens in 
a wide variety of FPT problems—either one of its first 
two moments can diverge. An abrupt phase transition 
is then observed— '^{k — F 0) = 00 and the introduction 
of unbinding is asserted to speed up turnover regardless 
of {Ton)- The asymptotic behavior of k'^jY high val¬ 
ues of {Ton) then depends on the tail of the catalysis 
time distribution {t — F 00), and it can be shown [49] 
that for /(t) ~ ^-(1+“) with 0 < a < 1 (1 < a < 2), 

C/T -- {Ton)~^ (C/T -- One example 

for this type of behavior is the above-mentioned case of 
diffusion mediated search for which a = 1/2. 

Universal behavior at fast binding times. When 
{Ton) approaches zero, “large fc” solutions to Eq. (3) are 
anticipated provided 'I'(A:) is asymptotically positive (see 
Fig. 2B). The behavior of 'I'(A;) in this limit is governed 
by the behavior of feat (t) at short times and we progress 
by assuming that feat (t) ^ (t ^ 1). Three differ¬ 
ent regimes, illustrated in Fig. 3, are then noteworthy 
[49] . When —1 < a < 0, 4'(A:) approaches zero from 
above as ^ k~^, and k^^Y ~ {Ton)~^- On the other 
hand, when a > 0, 'I'(A:) approaches zero from below as 
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~ —k~^, there are no “large /c” solutions to Eq. (3), and 
is finite (can be zero). 

The case a = 0 is a bit more delicate. Assuming feat {t) 
has a Taylor expansion near the origin we denote ujn = 


n\ 




and note that by construction 


= feat (0) > 0. Implementing a treatment similar to 
one given in the previous section we find for fc ^ 1 [49]: 


x=o/fc 


where i?oo = 


~ l + 2fl:oo(Xcx=/fc) + i^cx=(l + i^oo)(Xoo/fe)^ ’ 

= (w 2 - 2 a;^)/ 2 wi (compare with '^{k) 
for fc 1 above). A large k solution to Eq. (2) is then 
found only when fcatif^ decreases steeply enough near 
the origin, i.e., when f^at (0) < -flat (0) Xoo > 0, 
and is given by 


C/T - Xoo • (5) 

In particular, note that in this case kff^^ {Ton)~^^^- 
On the other hand, when < 0, we once again find 
that fc^/“|(ro„)=o is ^ constant. 

Conclusions. In this letter we took advantage of the 
Michaelis-Menten reaction scheme to provide a unified 
analysis of the optimal restart problem. The incorpora¬ 
tion of restart into an existing first passage time problem 
modifies its behavior. The mean first passage time then 
becomes a function of the restart rate and the question 
of optimality naturally arises. Here, we have developed a 
formalism which can be used in order to study the effect 
of restart on generic first passage time problems. A prime 
corollary of our study is the identification of two regimes 
at which the optimal restart rate displays universal be¬ 
havior in the sense that it is solely governed by a handful 
of key parameters (see Fig. 4 for illustration). The re¬ 
sults we have obtained are applicable to many fields. In 
particular, we note that randomized computer algorithms 
[54, 55[ often exhibit heavy tailed run time distributions 
[56, 57[ and restart could hence drastically improve per¬ 
formance in these cases and others [58, 59[. 

The optimal restart problem is intimately related with 
the idea of kinetic proofreading. Independently proposed 
by Hopfield [60[ and Ninio [29[, and studied by multiple 
authors since [61-69], the kinetic proofreading scheme 
suggests a way in which an enzyme can amplify small 
differences in the unbinding rates of two substrates—one 
right, the other one wrong, such that —in 

order to discriminate them with high fidelity. The basic 
idea is that by molding the catalysis time distribution, for 
example into something that resembles a sharp delay, the 
ratio between the right and wrong turnover rates can be 
made arbitrarily large. Consider now a case in which one 
is provided with a desired profile of the turnover rate as a 
function of the unbinding (restart) rate, and is then asked 
to conjure a catalysis (i.e., first passage) time distribu¬ 
tion that would yield this profile. For example, think of 



Figure 4: Color Online. Optimal unbinding (restart) 
rate vs. mean binding time for different catalysis time 
distributions. The black (full), blue (dashed), and red (dot¬ 
ted) curves correspond to kfff, and to asymptotic approxi¬ 
mations of kfff in the (To^) —>■ 4'o, and {Ton) 0, limits 
respectively. Distributions (A-Double Exponential, B-Triple 
Exponential, C-Gamma, D-Levy-Smirnov) are drawn at the 
bottom left corner of each panel [49]. Dashed blue lines in 
panels (A-C) are drawn according to Eq. (4). A phase tran¬ 
sition occurs in panel (D) where {Teat) = oo and asymptotics 
is governed by a ~ {Ton)~^ scaling law. Dotted red lines in 
panels A & B are drawn according to Eq. (5). A different 
behavior is found in panels C & D—in accord with the t —> 0 
asymptotics of fcat{t) and the scaling laws that are summa¬ 
rized in Fig. 4. Note how in panel (A) the exact solution 
coincides with the asymptotic approximations as expected. 


a scenario in which an enzyme wishes to select only the 
substrates whose unbinding rates fall within a small range 
(band pass filter), or above/below some cutoff (high/low 
pass filter). How could this be done? Quite surprisingly, 
a formal solution to this highly non-trivial problem can 
be readily obtained by solving Eq. (2) for fcat{koff) 
given kturn{koff) [fO]. This observation paves the way 
towards intelligent design of “Michaelian Filters”—a con¬ 
cept which we will further develop elsewhere. 
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1 Supplementary information for figures 

1.1 Details of distributions used in Fig. 2 

Turnover rates were calculated for catalysis time distributions with {Teat) = 1, 
and the mean binding time was taken to be (To„) = 0.1 in all cases. The 
following distributions were used (numbers match lines No. 1-6 in Figure 2): 

1. Log-Normal distribution: 


fcat{t) 


1 \ - {ln{t) - 

[ 2a2 



with a = \Jln{'i), and (3 = —ln{3)/2. The variance of this distribution is 

{Teat) =2. 

2. Double-Exponential distribution (a.k.a. Hyper-Exponential): 

fcat{t) = + (1 - , 

with fci = 0.5, fc 2 = 2, and p = 1/3. The variance of this distribution is 

G\Teat)=2. 

3. Weibull distribution: 

fcat{t) = aP {Pt)°‘~^ , 

with a = 0.7209, and /? = 1.232. The variance of this distribution is 

{Teat) = 2. 

4. Double-Erlang distribution (a.k.a. Hyper-Erlang): 

fcat{t) = pklte~'"^^ + (1 - p)klte~'"^* , 

with ki = 22.8781, k 2 = 1.4184, and p = 0.31. The variance of this 
distribution is a‘^{Teat) — 1.0612. 

5. Double-Exponential distribution (a.k.a. Hyper-Exponential): 

fcat{t) = pkie~'^^* + (1 - p)k 2 e~'‘^* 

with ki = 0.5, k 2 = 10/9, and p = 1/11. The variance of this distribution 
is a‘^{Teat) = 1-2 . 

6 . Exponential distribution: 

fcat{t) = fee”'"*, 

with k = 1. The variance of this distribution is (j'^iTeat) = 1 • 
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1.2 Details of distributions used in Fig. 4 

Maximal turnover rates were calculated for the following catalysis time distri¬ 
butions (by panel order): 

A. Double-Exponential distribution: 

fcat(t) = pkie~'"^* + (1 - , 

with k\ = 1, /c 2 = 10.25 and p = 0.1025. 

B. Triple-Exponential distribution: 

fcat(t) = pikie~^^* + p2k2e~'"^^ -b (1 - pi - P 2 ) , 

with fci = 1, /c 2 = 6, fca = 15, pi =0.1 and p 2 = 0.3. 

C. Gamma distribution: 

T(a) 

with a = 0.1942 and [5 = 1.0221. 

D. Levi-Smirnoff distribution: 

with D = 1 and t = 0.5. 
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2 Supplementary Fig. SI 



Asymptotics of the optimal restart problem, at {Ton) ~ is gov¬ 
erned by the first three moments of the catalysis time distribu¬ 
tion. Distributions No. 1-4 are all of type Double Erlang with (Teat) = 1 
and (j{Toat) = "v/S- They differ, however, by having increasing skewnesses— 
i?o = {0.7, 1.5, 3.5, 14.8} respectively (see below). Inset. When (Ton) ~ Tq, 
kmax is linear in A = {To„) ^ and inversely proportional to Rq. As Rq 

increases, the maximizing unbinding rate becomes lower as to avoid premature 
termination of the catalytic process. 

2.1 Details of distributions used in Fig. SI 

Distributions No. 1-4 are all of type Double-Erlang 

fcat{t) = pklte~^^* + (1 - p)klte~'^^*. 

The following parameters were used: 



P 

ki 

k2 

(Teat) 


i?0 

1 

1.93* 10-1 

0.55 

5.437 

1 

4 

0.670 

2 

1.23* 10-1 

0.45 

3.875 

1 

4 

1.467 

3 

4.94* 10-^ 

0.3 

2.833 

1 

4 

3.477 

4 

4.6*10-^ 

0.1 

2.19 

1 

4 

14.76 
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3 Supplementary information for text 

3.1 Derivation of Eq. 3 in the main text 

Before we derive Eq. 3, and the statements which follow it, let us hrst obtain a 
preliminary result. Unbinding is said to be excitatory whenever dktum/dkoff > 
0. Utilizing the dehnition of ktum (Eq. 2 in the main text), we see that this 
condition is equivalent to 


fcat{k,ff ) ( 1 - /c.t(fco//)) + fco// (1 + fco// (T„„)) 

{{Ton) koff + 1 — fcat{koff ) j 


> 0 . 


(3.1) 


Noting that the denominator on the left hand side of Eq. (3.1) is always positive, 
we multiply both sides by it and rearrange to get 


Ltikoff) (l - fcatikoff)) > -koff (1 + koff (Ton)) • (3-2) 

The Laplace transform of a probability density function is monotonically 
decreasing in koff. We can thus take Eq. (3.2) and divide both sides by 
-dfcat{koff)/dkoff to get 


fcat{koff) (^1 fcat{koff 

d Tat {koff) 
dkoff 


^ koff (1 + koff (Ton)) ■ 


(3.3) 


Further noting that the unbinding rate is always positive, we rearrange Eq. (3.3) 
to show that the following condition 


Jcaty^of f ) (jcaty^of f) 

dfoat{k„ff) ,2 

dkoff off 


koff 


> (To 


is equivalent to the existence of excitatory unbinding. 
Examining the left hand side of Eq. (3.4), we dehne 


and conclude that 


l'(fc) = 


fcat{k) {jcat(k) - 1 


dfoatik) 

dk 


/c2 


1 


(3.4) 


(3.5) 


^{koff)>(Ton)^'j^>(^. (3.6) 

dkoff 

A similar set of steps can then be applied to show that in the case of inhibitory 
unbinding one gets 

nkoff)<(Ton)^^<Q, (3.7) 

dkoff 

as illustrated in the hgure below. 
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k. 


turn 



1.05 * ^^tum 


^^tum 


Excitatory 

unbinding 


Inhibitory 

unbinding 


V(fco//^)><r„„> 'i’ikoff) < {Ton) 


0 . 8 , 


0 


0.5 


1 I, 1.5 2 2.5 3 



max 


Unbinding rate 


Replacing inequalities with equalities in the above derivation, one can readily 
convince himself that 



(3.8) 


i.e., the stationary points of ktum are the solutions of 'ii{koff) = {Ton }—as is 
stated by Eq. 3 in the main text. We furthermore note that a stationary point 
of kturn for which < 0 {'^'{k) > 0) is a a local maximum (minimum) 

of the turnover rate. Indeed, the amalgamation of Eqs. (3.6) and (3.7) asserts 
that dkturn/dkoff is positive (negative) if and only if ^'(fcojj) — (Ton) is positive 
(negative). Since at a stationary point ^'(A:) — {Ton) = 0, 4>'{k) < 0 {'^'{k) > 0) 
implies that at this point dktum/dkof f transitions from being positive (negative) 
to being negative (positive), thus resulting in a local maximum (local minimum). 
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3.2 Approximating 

3.2.1 A general recipe 

Before obtaining approximations for at particular limits, we give a recipe 
illustrating how this could be done in general. Consider the limit k —> kum- 
An approximation for at this limit, can be found by taking the following 

steps; 

1. Find the third order power expansion of the Laplace transform fcat{k) in 
the limit k —> ku^ (going to third order is required in order to capture 
the hrst order behavior of 5'(fc)). If the Laplace transform is not analytic 
move to step number (5) below. 

2. Replace the third order power expansion of fcat(k) with the following 
rational function approximation 


fcat{k) 


1 + ak 
1 + bk + ck'^ 


(3.9) 


Choose the parameters a, b, and c such that the third order power expan¬ 
sion of the right hand side of Eq. (3.9) agrees with that found in step (1) 
above. 


3. Substitute fcat{k) in the dehnition of ^'(/c) (see Eq. 3.5 above) with its 
approximation as given by Eq. (3.9). 

4. Note that, under the approximation made in step (3), ^'(/c) = (To„) (Eq. 
3 in the main text) is a quadratic equation. Solve it. 

5. In cases where a power expansion of fcat{k) does not exist, leading order 
behavior must be extracted in a different manner. In some cases, approx¬ 
imate solutions to Eq. 3 in the main text can still be found. Examples 
are given below. 


3.2.2 Approximating 4t(/c) in the fc —> 0 limit 

1. Power expansion. When all moments of the catalysis time distribution 
are finite fcat{k) can be written in the following form [1] 

fcat{k) = Y,^{-kr , (3.10) 

where M„ = feat is the moment of the distribution. Truncating 
the infinite series in Eq. (3.10) to obtain a third order power expansion we have 

fcatik) ~l-Mik + ^k^-^kf^ . (3.11) 

Z D 


2. Rational function approximation. Note that to third order 

l+gfc ^ 
l+6/c+cfc2 

1 — (& — a)k -\- (b^ — ab — c)k^ — {b^ — ab^ + ac — 2bc)k^ . 


(3.12) 
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Equating corresponding coefficients for the power expansions that appear on 
the right hand side of Eqs. (S.ff) and (3.12), one can write equations for the 
set {a, b, c} and solve them to obtain 

_ 6Mi(M2 - Ml) - M3 3M1M2 - M3 _ 3 M| - 3M1M3 
“ 3(2M2 - M 2 ) ’ 3(2M2 - Ma) ’ ^ 6{2Mf -M 2 ) ' ^ ^ 

Substituting these parameters into the right hand side of Eq. (3.9) we obtain 
the following rational function approximation 


fcatik) ~ 


3{2Mf-M2) j 


1 + 


{ 3M\ M 2 






\3(2Mf-M2)'} ^ + ( 6 ( 2 M 2 -M 2 ) J 
3. Approximating 5'(fc). Substituting Eq. (3.14) into Eq. (3.5) we get 


(3.14) 


1 

'i'(fc) 


2 Mi 


M2-2M{ 


+ 2 


2M1M3-3MI 

3[M2-2Mfj^ 


(6iV/(^ — (iM2 M) (2iV/i Ad 2 — 3 Ad 2 

9(M2-2Mf)3 


k 


k\ 


(3.15) 


Defining 



_M2-2M2 

- 2Mi ’ 

(3.16) 

and 

_ 2 M 1 M 3 - 3M| 

" ~ 3 (M 2 - 2M2)^ ’ 

(3.17) 

we rewrite Eq. (3.15) in 

a more compact form 


^(fc) ^ 4-0 ( 

'1 + 2 i?o i^ok) + [i?o(l + ^ 0 )] (^-ofe)") . 

(3.18) 


Examining Eq. (3.18), it can be seen that its right hand side depends on 
two parameters only 4^0 and Rq. We further note that ^'(fc = 0) = and that 
^'o defines a time scale by which we can normalize other times in our problem. 
Defining 


t = fc = ^<0 «'(fc) = •^o'^normiK) , (3.19) 

we can write a normalized version of Eq. (3.18) 

^'norm('^) — (l + + [^o(l + ^o)] ) (3.20) 

in which was eliminated, and Rq —the only remaining statistic—is dimen¬ 
sionless. 


4. Solving the governing equation. Substituting 4'(fc) in Eq. 3 (main 
text) with its approximation, as given by Eq. (3.18), we have 



A. 


B. 


■ fli 




A 

-1 < Ao„ < 0 

1 

1 

Ro > 0 

1 

1 

1 

1 

-Kn 

1 

1 — rJ --- - > 



1 + 2_Ro (^'ofc) + [-Ro(l + -Ro)] — ^ 0 / (Ton) ■ (3-21) 

Equation (3.21) is based on a small k approximation of and we will hence 
try to solve in this regime. Small k solutions to Eq. (3.21) can only be found in 
the (Ton) —>■ 5'o limit and we therefore see that, since (Ton) > 0 by dehnition, a 
requirement for their existence is 5'o ^ 0 • Two different scenarios should then 
be distinguished based on the hrst order (linear) behavior of the left hand side 
of Eq. (3.21). If i?o > 0, a small k solution exists only for 0 < (Ton) < Tq. On 
the other hand, if i?o < 0, a small k solution exists only for (Ton) > ^'o- The 
two solutions of Eq. (3.21) are 


k 


± _ 
ext 


1 

(1 + Ro) ^'o 



To 

(Ton) 



(3.22) 


but > 0 in the limit (Ton) —t Tq, while does not. We therefore choose 

k'^^-i- and continue our analysis. 

Utilizing Eqs. (3.19) and (3.22) we write k'^^^ in a normalized form 


K 


+ 

ext 


1 


(1 + Ro) 


Ro — An 


1 + An 


Rq 


- 1 


(3.23) 


where we have defined Ao„ = {(Ton) — To) /Tq. Recalling Eq. (3.20), noting 
that 0 < (Ton) < To — 1 < Aon < 0, and that (Ton) > Tq Aon > 0, we 
schematically illustrate the difference between the i?o > 0 case and the 7?o < 0 
case in the hgure above. Most importantly, we note that from our foldings in 
section (3.1) it follows that when Ro > 0 {Rq < 0), is a local maximum 
(minimum) of the turnover rate. 


5. Special cases. The analysis given above was based on a power expansion 
of fcat{k) in which we assumed that the first three moments of the catalysis 
time distribution are hnite. In particular, 5'o was dehned based on the hrst two 
moments of the distribution and is hence not well dehned when these diverge. 
Divergence of the hrst two moments can occur if the PDF has a “heavy tail” as 
happens, for example, when 

fcat{t) — ^ 1 +a ’ (3.24) 
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for i 2> 1 and 0 < a < 2. However, even when this is the case, one can still say 
something about the behavior of 5'(fc) at small k, and we distinguish between 
two different cases. 

When 0 < a < 1, the first and second moments of the catalysis time distri¬ 
bution both diverge. A small k approximation of the Laplace transform is then 
given by [1] 

f,at{k) - 1 - , (3.25) 

a 

where r(a;) is the Gamma function. Substituting Eq. (3.25) into Eq. (3.5) we 
see that to first order in this limit 


— -1. (3.26) 

a k 

Neglecting lower order corrections, we see that 5''(fc) < 0 and conclude that as 
{Ton) oo the following asymptotics holds 


k 


max 

off 


I-a 1 
a {Ton) ' 


(3.27) 


When 1 < a < 2, the first moment of the catalysis time distribution is 
finite but the second moment diverges. A small k approximation of the Laplace 
transform is then given by [1] 

foat{k) - 1 - {Toat) k + . (3.28) 

a(a — 1) 

Substituting Eq. (3.28) into Eq. (3.5) we see that in this limit 


^'(fc) 


C7-r(2-a) 1 

a {Toat) ^ 


(3.29) 


Neglecting lower order corrections,we see that 'E''(/c) < 0 and conclude that as 
{Ton) OO the following asymptotics holds 


[ C-r{2-a) 

\a {Toat) {Ton)] 


(3.30) 


3.2.3 Approximating 'I'(/c) in the k ^ oo limit 

1. Power expansion. When the PDF of the catalysis time distribution is 
analytic in the vicinity of t = 0 we can Taylor expand it 


fcat{t) = Y] 

n=0 



n\ 


oo 

y ^ ^n+1 

n=0 


(- 0 " 

n\{n -I- 1)! ’ 


where we have defined 


(3.31) 




= n\ 



from reasons that will soon become clear, 
series can then be written as 


fit) 


(3.32) 


t=oJ 

The Laplace transform of such a 
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Truncating the infinite series in Eq. (3.33) to obtain a third order power expan¬ 
sion we have 

t ^^2 1^3 1 /Q Q .X 

(3.34) 

2. Rational function approximation. Note that to third order 


1 + ak a / 1 

1 -|- &fc -|- cfc2 c \ k 


a/l\ ab — c / 1\ ^ ab'^ — ac — be / I 


Equating corresponding coefficients for the power expansions that appear on 
the right hand side of Eqs. (3.34) and (3.35), one can write equations for the 
set {a, 5, c} and solve them to obtain 

QuJi (uJ2 - 2 ujI) 2uJ3-Q(^1UJ2 6{uj2-2ujI) 

a — ——^ ^ ^ , 0 — - —c — —^^ . (3.36) 

2idiOJ3 — 3cJ2 2(jJiOJ3 — ^^2 2uJiOJ3 — 3cj2 

Substituting these parameters into the right hand side of Eq. (3.9) we obtain 
the following rational function approximation 


fcat{k) ~ 


^ k Gu^-2^) \k 

\ 2cc;iCc;3—3cj2 / 

1 I f u : f ^(^^- 2 ^ 1 ) \ u2 

y — J I 2aJia;3 —/ 


3. Approximating ^(/c). Substituting Eq. (3.37) into Eq. (3.5) we get 


1 2 a;i 7,2 i o 2coiu)3—3u)2 

^{k) - [u;2-2u;^\ ^ ^ ^ ^{u^2-2u^lY 


Defining 


(6cc>^ — 6 co>2^i ) (2cJia;3 —) 

9(cd2—2cd^)3 


_ OJ 2 — 

Xoo = ? 

2uji 


_ 2 cjiW 3 — 3 w| 

3(cU2- 2 ^ 2)2 

we rewrite Eq. (3.38) in a more compact form 


'l'(fc) - {{xjkf + 2 Roo (Xoo^fc) + ^oo (^oo + 1 )) 
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B. 





'viKff) 



Examining Eq. (3.41), it can be seen that its right hand side depends on 
two parameters only Xoo and Roo- We further note that in the limit k —> oo, 
Xxk~'^, and that x^ dehnes a time scale by which we can normalize 
other times in our problem. Dehning 

t = k = XqK, 4'(A:) = Xq^'^N orra{K) (3.42) 

we can write a normalized version of Eq. (3.41) 

^ N orm {K)~{K^ + 2R^K + R^{R^ + l))\ (3.43) 

in which Xoo was eliminated, and f?oo—the only remaining statistic—is dimen¬ 
sionless. 

4. Solving the governing equation. Substituting 4'(fc) in Eq. 3 (main 
text) with its approximation, as given by Eq. (3.41), we have 


_Xcxj/fc^_ 

1 -I- 2R^Xoo/k + Roo [Roo + 1) {Xoo/ky 


(Ton) ■ 


(3.44) 


Equation (3.44) is based on a large k approximation of 4'(fc) and we will hence 
try to solve it in this regime. Large k solutions to Eq. (3.44) can only be found 
in the {Ton) 0 limit. Two different scenarios, schematically illustrated in the 
figure above, should then be distinguished based on the leading order behavior 
of the left hand side of Eq. (3.44). 

When Xoo > 0, lim 4'(fc) = 0+, i.e., 4'(/s) > 0 as it approaches zero from 

k—^oo 

above. The two solutions of Eq. (3.44) are then 


^ext Xo 


Xoo (To] 


-R oo -R O 


(3.45) 


but —>■ oo in the limit (To„) —> 0, while k~^f does not. We therefore choose 

and continue by noting that since 4''(A:) < 0 in this limit the Endings in 
section (3.1) assert that k^^^ is a local maximum of the turnover rate. Utilizing 
Eqs. (3.42) and (3.45) k^^^ can also be written in normalized form 


Kt.t 



(3.46) 
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where we have dehned Son = Xca {Ton) ■ 

When Xoo < 0, lim = 0“, i.e., <0 as it approaches zero from be- 

k—^oo 

low. Since (Ton) > 0 by dehnition, there are no large k solutions to Eq. 3.44 in 
this case. It is however possible, that 5'(fc) is positive for smaller values of k, in 
which case, the continuity of and the fact that it is negative at large k guar¬ 
anties that there is a hnite value of k, kt^nesh, for which ^(kthresh) = 0. This 
translates (as visualized in the figure above) to kmax {{Ton) —> 0) —> kthresh ■ 


5. Special cases. The analysis given above was based on an expansion of 
fcat{k) in powers of k~^. In particular, when defining Xoo we have assumed 
that 0 < fcat{ 0 ) = wi < oo. In contrast, let us now consider a power law 
behavior of the PDF 

foat{t)~toir, (3.47) 

where t ^ 1, wi > 0, a ^ 0, and a > —1 (the PDF cannot be normalized for 
a < —1). Noting that the large k asymptotics of the Laplace transform is then 
given by 

nOO 

foat{k)~ / wit“e-'^*df = wir(l + a)fc-(i+fo, (3.48) 

Jo 

we substitute Eq. (3.48) into Eq. (3.5) and find 


T(fc) 


—a 1 
a -h 1 fc ’ 


(3.49) 


to first order in this limit. 

Examining Eq. (3.49), we see that two different cases should be told apart. 
When — 1 < a < 0, /cat(t) diverges in the limit t —> 0, T(fc) > 0, and T'(fc) < 0. 
Neglecting lower order corrections, we conclude that as (Ton) —> 0 the following 
asymptotics holds 


•Jf ^ a+1 (Ton) ■ 


(3.50) 


On the other hand, when a > 0, fcat{t) zeros in the limit t —> 0, and 
5'(fc) < 0. Since (Ton) > 0 by definition, there are no large k solutions to Eq. 
3.44 in this case. It is however possible, that 4t(/c) is positive for smaller values 
of k, in which case, the continuity of 5'(fc) and the fact that it is negative at large 
k guaranties that there is a finite value of k, kthresh, for which ^ (kthresh) = 0. 
In this case, we once again find k'^j^ {{Ton) —> 0) —> kthresh ■ 
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